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We study the pair-breaking effect of rotation on a cold Fermi gas in the BCS-BEC crossover region. 
In the framework of BCS theory, which is supposed to be qualitatively correct at zero temperature, 
we find that in a trap rotating around a symmetry axis, three regions have to be distinguished: (A) 
a region near the rotational axis where the superfiuid stays at rest and where no pairs are broken, 
(B) a region where the pairs are progressively broken with increasing distance from the rotational 
axis, resulting in an increasing rotational current, and (C) a normal-fluid region where all pairs are 
broken and which rotates like a rigid body. Due to region B, density and current do not exhibit any 
discontinuities. 
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The surprising properties of superfiuids become most 
evident if one looks at rotating systems. But the rotation 
does not only reveal the superfluidity, it can also destroy 
it. To give an example, in nuclear physics, the strong 
reduction of the nuclear moment of inertia compared to 
its rigid-body value is a direct consequence of superflu- 
idity due to pairing correlations. But with increasing 
angular momentum, the pairing correlations are progres- 
sively destroyed and the moment of inertia increases to 
its rigid-body value. This pair-breaking effect of rotation 
was studied many years ago 

In trapped atomic Fermi gases, the picture is somewhat 
different, since, contrary to the situation in atomic nuclei, 
the coherence length is much smaller than the system 
size. It is therefore possible to create quantized vortices 
or even vortex lattices 0] , which allow the system to stay 
superfiuid while rotating. 

However, in a recent paper by Bausmerth, Recati, and 
Stringari Q it has been argued that it may be possible 
to put a trapped Fermi gas adiabatically into rotation 
without creating vortices. In that paper, the destruction 
of superfluidity by rotation is described in a way which 
is very different from the nuclear physics case: Instead 
of decreasing the value of the pairing gap with increas- 
ing angular velocity, the authors assume that the system 
separates into a paired and an unpaired phase, while the 
properties of the paired phase itself are not affected by 
the rotation. The authors consider the unitary limit, 
where the energy densities of the paired and unpaired 
phases are known from Quantum-Monte-Carlo (QMC) 
simulations . The phase boundary between the paired 
and the unpaired phases is determined by energy mini- 
mization: Near the rotational axis, the system prefers to 
stay superfiuid, i.e., to stay at rest, since the paired phase 
has a lower energy density than the unpaired one. But 
beyond a certain distance from the rotational axis, the 
centrifugal energy which the system could win if it partic- 
ipated in the rotation becomes equal to the energy which 
is needed to break the pairs. Hence, the non-rotating 
superfiuid core is surrounded by a rotating normal-fluid 
phase. At the interface separating the two phases, the 



density and the current are discontinuous. 

This picture is very intuitive, but it is lacking the mi- 
croscopic understanding of the pair-breaking mechanism. 
In the present paper we will therefore describe the ro- 
tating Fermi gas in the framework of BCS theory. The 
rotation is most easily described in the rotating frame, 
where the hamiltonian H (minus the chemical potential 
^ times the particle number N) is given by 
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+ 5^|(r)V3|(r)V3i(r)V3T(r)] , (1) 

where ^ is the Fermion field operator with components 
for (pseudo-)spin up (|) and down (j), m is the atom 
mass, p = — i?iV and L = r x p are momentum and angu- 
lar momentum, respectively, V{r) = m{uj'^z'^ + uj'^r'^)/2 
is the axially symmetric trap potential and 5 < is the 
coupling constant. The system is supposed to rotate with 
angular velocity Q around the symmetry (z) axis of the 
potential. 

If the system is large enough, such that the coher- 
ence length is small compared with the oscillator length 
associated with the trap potential, we can make use 
of the Thomas-Fermi (TF) or local-density approxima- 
tion (LDA), which amounts to treating the system at 
each point r as uniform with a local chemical potential 
A'ioc(r) = A* ^ ^(i")- Then p becomes a number instead 
of an operator, and the "cranking" term flL^ can conve- 
niently be written as fiL^ — v(r) • p, where v(r) = O x r 
is the velocity field corresponding to a rigid rotation. All 
quantities depend only parametrically on r via yLt;oc(r) 
and v(r). 

The gap, density, and current can all be derived from 
the normal and anomalous Matsubara Green's functions 
Q and JFt 5*1 . They have to satisfy the Gorkov equations, 
which in the presence of the cranking term QLz become 

{ihujr, -C + ^r■p)g + AT^ = h, (2) 

{ihujn+^ + v -P^J"^ + = 0. (3) 
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where we introduced the abbreviation ^ = ^(r, p) = 
p^/ (2m) — ^ioc(j), ujn denotes a fermionic Matsubara fre- 
quency, and A(r) is the gap. Note that we are neglecting 
the Hartree mean field, but anyway it would not qual- 
itatively change our results in the BCS-BEC crossover 
regime 0- Eqs. ([2]) and ([3]) can readily be solved for 
Q and J-. They are formally similar to those describing 
pairing between particles with unbalanced populations 
(see, e.g., Ref. Q), except that here the chemical poten- 
tials for the two spins are equal and the asymmetry is 
between states with opposite momenta (p and — p). 

In the case of a system without superfluid flow (like in 
our axially symmetric trap, as long as there are no vor- 
tices), the gap can be assumed to be real (A = A*). The 
gap equation is obtained in the usual way by summing 
T over ojn and integrating over p, with the result 

(4) 

where we defined the quasiparticle energies E± = E ± 
p • V, with E = v^2-fA2, and f{E) = I/^c^/C^bT) _^ 
denotes the Fermi function, T being the temperature and 
ks the Boltzmann constant. In Eq. ([3]), the divergence of 
the gap equation due to the contact interaction has been 
regularized in the usual way by expressing the coupling 
constant g in terms of the s-wave scattering length a Q . 

We are mainly interested in the BCS-BEC crossover 
regime, where it is known that the BCS description fails 
at higher temperatures, and in particular the BCS pre- 
diction for the critical temperature Tc is much too high. 
However, at zero temperature, BCS theory gives a rea- 
sonable description throughout the crossover. We will 
therefore restrict ourselves to the zero-temperature case, 
in which the Fermi function reduces to a step function, 
f{E) = e{-E). Hence, the factor [1 - f{E+) - f{E^)] 
is equal to 1 if both E+ and E- are positive and oth- 
erwise (at most one of the two energies E+ and E- can 
be negative). In other words, states with E± < are 
excluded from pairing. In order to better understand the 
role of these states, let us look at the occupation numbers 
p(r,p), which are obtained by summing Q over 

p{r, p) = i (l - I) [1 - fiE+)] + ^ (l + I) f{E-) . (5) 

For states with both E+ > and E^ > 0, this reduces to 
the usual BCS expression. But if a state with momentum 
p has E^ < 0, its occupation number is equal to 1. The 
corresponding time-reversed state with momentum — p 
has then E^ < and its occupation number is equal to 
0. As we will see below, this gives rise to a normal-fluid 
(rotational) current. 

It is easy to see that the energies E± can only become 
negative if the velocity v exceeds a critical value such 
that 

p'fV> A. (6) 

Here we have introduced the abbreviation p'p = 
^/2rnJ/^, where fi^^^ = fiioc + mv'^/2 denotes the local 



chemical potential which includes the effect of the cen- 
trifugal force, and p'p is the corresponding local Fermi 
momentum. For a given z coordinate, the condition ([5]) 
is fulfilled beyond a certain distance r±i{z) from the ro- 
tational axis, since the velocity increases as v — Qr±. At 
smaller distances, the energies E± are always positive, 
i.e., the system is in the usual superfluid phase and does 
not participate in the rotation. Beyond rj_i, the gap is 
reduced by the rotation. We will call this region, where 
a rotational current exists although the gap is non-zero, 
the partially paired phase. Finally, at a certain distance 
r±2, the gap vanishes and the system enters the normal 
phase where it rotates like a rigid body. 

If the condition ^ is fulfilled, i.e., for 7-^ > rj_i, one 
can easily see that the energies E± can become negative if 
the momentum lies between two limits p_ and p+ which 
are given by 

pI = p'^ + ± 2m^p'^v^ - , (7) 

The integrand of the gap equation Q is only affected by 
the rotation if p lies between p_ and p+ . Integrating Eq. 
Q over the angle between p and v and dividing both 
sides of the equation by A, we obtain 

-^[r*(^-H-r*(s-9i<«) 

The first integral is the same as in the gap equation with- 
out rotation while the second one is the contribution of 
the f{E±) terms due to the rotation. 

In the weak-coupling limit, when A <C /ijoc, the pair- 
breaking effects appear already at extremely low angular 
velocities f2. In this case it is possible to evaluate the in- 
tegrals in Eq. ([8]) analytically, and one can show that the 
critical velocity for which the gap disappears is given by 
Vc — (e/2) Ai,=o/pf, where e — 2.71 . . . denotes Euler's 
number. Hence, for a given z coordinate, the radial co- 
ordinates r_i_i^2 separating the fully paired from the par- 
tially paired and the partially paired from the unpaired 
phase, respectively, are the solutions of the equations 

PF[r^i,z)VLr^i = An=n[r±_i,z) , (9) 

PF(f_L2, z)Q.r±_2 = -j^Ai-i^o{r±_2, z) . (10) 

In the crossover regime, the situation is more compli- 
cated, since the gap A may be comparable with ^loc- 
Therefore the integrals have to be evaluated numerically. 
In addition, the rotation can now be much faster and the 
centrifugal force can lead to a sizeable change of the den- 
sity profile and it is necessary to readjust the global chem- 
ical potential /i as a function of in order to keep the to- 
tal number of particles fixed. The density per spin state, 
p(r), is obtained by integrating the occupation numbers 
over p. Using Eq. ([5]), one obtains 

(11) 
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FIG. 1: From top to bottom: density per spin state p, gap 
A, and current |j| in a rotating Fermi gas (4 • 10"" atoms in 
an isotropic trap with frequency to) in the BCS phase as a 
function of the distance r± from the z axis, for z — 0. The 
solid lines correspond to a gas rotating with angular velocity 
il = 0.05u>. For comparison, the results for the non-rotating 
case (ground state) are shown as the dashed lines. 



The second term arises from the f{E±) terms and ex- 
ists only if the condition ([6]) is fulfiUed, i.e., beyond 
r±i. Between r±i and r±2, the density goes smoothly 
from its value with pairing to the value without pairing, 
limA^o p(r) — p'p{v)/{QTi'^'h?). Once we have calculated 
the density, we can obtain the total number of particles 
by integrating the density over space. This allows us to 
determine the value of the chemical potential. 

An interesting quantity is the current density, which 
can be obtained by multiplying the occupation numbers 
with p/m and integrating over p. From Eq. ([5]) it is clear 
that for < i, i.e., close to the rotational axis where 
the condition ([6]) is not satisfied, the current vanishes as 
it should in the superfluid phase. Beyond r^i, the result 
can be given in closed form as 



67r2n3 



2\3/2 



(12) 



One sees that in the partially paired phase the current 
increases with decreasing gap and it correctly approaches 
its rigid-body limit if one approaches the unpaired phase: 
limA^oj(r) = /o(r)v(r). 

Let us now discuss some numerical results. We con- 
sider a system with iV = 4 • 10^ atoms (2 • 10^ atoms per 
spin state) in two cases: (a) close to the BCS limit, with 
l/{kpa) — —2 [kp = Pf{y = 0)//i], and (b) at unitarity, 
i.e., in the limit a — > cxd. We do not consider the BEG side 
of the cross-over, since as soon as the chemical potential 



FIG. 2: Same as Fig. [T] but for a unitary Fermi gas rotating 
with angular velocity = 0.45cj. 



becomes negative, the energies E± are always positive, 
i.e., the molecules in the BEG phase are never broken 
by the rotation. For simplicity we choose a spherically 
symmetric trap {ujz — wj_), but this will not qualitatively 
change our results. In the figures, we will use the har- 
monic oscillator units set by the trap potential, i.e., tnj 
for energies and Iho — \J hjiraiS) for lengths. 

Let us first discuss the BGS case. In this case the pair- 
ing is so weak that it does not appreciably influence the 
density (upper panel of Fig.[T|). It is also very fragile, i.e., 
the moment of inerticL which can be calculated within lin- 
ear response theory [S], must be measured at extremely 
low angular velocity. Already for an angular velocity as 
small as ri = O.OSw, the gap (second panel of Fig. [1]) is 
zero in a large part of the system. Because of the small 
angular velocity, the centrifugal force has no effect on the 
density, either. Looking at the gap, one can clearly see 
the point r±i{z = 0) = 5.1 Iho where the results for the 
non-rotating (dashed line) and the rotating (solid line) 
system start to differ, and the point r±2 =6.1 Iho where 
the gap goes to zero. The three regions are even more 
evident in the current (lower panel of Fig. [l}: The cur- 
rent starts to be non- vanishing at rj_i and it has a kink 
at r±2 where it reaches the rigid-body value. 

More interesting are the results in the cross-over 
regime, where the gap is strong enough to support a rela- 
tively fast rotation. In Fig. [2] we display the density, gap 
and current (from top to bottom) for a system at the uni- 
tary limit rotating with = 0.45 (solid lines; for com- 
parison, the density and gap of the corresponding non- 
rotating system are shown as the dashed lines). In this 
case, the centrifugal force leads to an oblate deformation 
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of the system: The chemical potential /i and the axial s ize 
of the system, which is determined by z^ax = •\/2/x/m/a;, 
decrease (in the present example, n decreases from 81.7 
to 81.2 huj), while the radial size, which is determined by 
r±max — •\/2/i/[m(cj^ — ^^)], increases. The increase of 
the radial size is visible in the upper panel of Fig. [21 where 
the density is shown as a function of r± for z = 0. The 
depletion of the density in the center is a consequence of 
the reduced chemical potential. This is also the reason 
why the gap in the center decreases with the rotation. 
Due to the strong pairing, the gap has a direct effect on 
the density. This is the reason for the kink in the density 
profile at r± = r±2- However, we stress that the density 
stays continuous at r±^2- 

The fact that, in contrast to the results of Bausmerth 
et al. 3, the density, the gap, and the current remain 
continuous functions of r±_ is the main statement of the 
present paper. In fact, if we followed the arguments given 
in Ref. we would find a similar discontinuity as they 
do. The only difference with their result would be the 
different numerical values of the parameters and ^jv 
which determine the relationship between the density and 
the local chemical potential [jiioc = ^sfi.^(67r^p)^/'^/(2m) 
or ^'i^^ = ^Ar?i^(67r^p)^/^/(2m) in the superfluid and nor- 
mal phase, respectively]. In BCS theory without mean- 
field shift, one obtains = 0.59 and £,n — 1, whereas the 
QMC resuhs used in Ref. [1] are £,s = 0.44 and £,n = 056. 
If one excluded the possibility of an intermediate "par- 
tially paired" phase, as in Ref. [sj, the system would 
have to split into a fully paired superfluid and a fully 
unpaired normal-fiuid phase, and the density would have 
a discontinuity across the phase boundary with a ratio 
Pn/ps = fe/^Ar)^/^ which gives 0.73 with the BCS 
results and 0.85 with the QMC results for and ^n- 
From this we see that, even if BCS theory is not capa- 
ble to give the right numbers for and ^at, the ratio is 
semi-quantitatively correct. Anyway, even if our results 
for the unitary limit might not be very precise, we be- 
lieve that they are qualitatively correct and that between 
the ordinary normal and superfluid phases there will be 
a region in which some pairs are broken while others stay 



unbroken. In particular, we checked that in the region 
between rj^i and r±2 our energy density is lower than 
both that of the non-rotating superfluid and the ridigly 
rotating unpaired gas. 

We emphasize that the existence of the intermediate 
region is not a finite-size effect, but it survives in ar- 
bitrarily large systems. For instance, if the trap was 
a flat potential well instead of a harmonic oscillator, 
the ratio of the two radii r±2 and r^i would become 
r±2/r±i = e/2 = 1.36 according to Eqs. ^ and (fTOl) . 
independently of the size of the system and of the an- 
gular velocity of the rotation. In the harmonic oscillator 
the intermediate region is smaller since the gap decreases 
with increasing r± already in the non-rotating case. 

This does not mean that finite-size effects do not play 
any role. For instance, the abrupt decrease of A for 
r_L rj_2 is an artefact of the TF approximation, which 
requires that all spatial variations be slow compared with 
the length scale set by the coherence length. A necessary 
condition for this is A ^ huj. In a true quantum cal- 
culation, the profiles of A, p and |j| would be rounded 
and no sharp interface between the different phases could 
be defined. In addition, beyond a certain critical angular 
velocity ilc the gap should completely disappear, even on 
the rotational axis (Toj . 

An interesting extension of the present work is to study 
a system which is deformed in the xy plane, i.e., in the 
plane perpendicular to the rotational axis. This question 
is very important since it is impossible to put the sys- 
tem into rotation without such a deformation (of course, 
once the system rotates, the deformation can be switched 
off and the conservation of angular momentum ensures 
that the system keeps rotating). In the deformed case, 
also the superfluid part of the system has a non- vanishing 
current, with an irrotational velocity field. Another im- 
portant question concerns the collective excitations of 
the rotating system, in particular the radial quadrupole 
mode whose precession is used in current experiments for 
measuring the angular momentum of the system [ll| . In 
order to stay in contact with the experiments, tempera- 
ture effects should be taken into account, too. 
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